Entanglement requirements for implementing bipartite unitary operations 
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We prove, using a new method based on map-state duality, lower bounds on entanglement re- 
sources needed to deterministically implement a bipartite unitary using separable (SEP) operations, 
which include LOCC (local operations and classical communication) as a particular case. It is known 
that the Schmidt rank of an entangled pure state resource cannot be less than the Schmidt rank 
of the unitary. We prove that if these ranks are equal the resource must be uniformly (maximally) 
entangled: equal nonzero Schmidt coefficients. Ifigher rank resources can have less entanglement: 
we have found numerical examples of Schmidt rank 2 unitaries which can be deterministically im- 
plemented, by either SEP or LOCC, using an entangled resource of two qutrits with less than one 
ebit of entanglement. 

PACS numbers: 03.67.Ac 



I. INTRODUCTION 

It is possible to carry out nonlocal quantum op- 
erations on multipartite systems using only local 
quantum operations and classical communications 
(LOCC) provided that the parties involved have ac- 
cess to a suitable entangled state, referred to as a 
resource. Given a large enough resource it is always 
possible to use teleportation to send all inputs to 
one party, who performs the operation and then dis- 
tributes the results to the other parties using tele- 
portation. In some cases it is possible to perform 
a nonlocal operation with less entanglement than is 
required by teleportation [iHEl- The question then 
arises as to how much entanglement is really neces- 
sary in order to implement a given nonlocal opera- 
tion. 

Our first result, that the Schmidt rank of the re- 
source must be at least as great as that of the unitary 
[Theorem 1(a)], follows rather immediately from the 
fact that it is a separable (SEP) operation. This is 
analogous to the result given in [3] in which prob- 
abilistic (i.e. SLOCC) implementations are consid- 
ered. Since SEP is contained in SLOCC, our The- 
orem 1(a) can be seen as a consequence of the re- 
sult in [7| , however we provide an independent proof 
along the way to our main result. 

In contrast to the probabilistic case, the determin- 
istic implementation of a unitary is only possible if 
the state meets certain entanglement requirements. 
For one thing, the entanglement of the resource 
must be at least as great as the entangling power 
of the unitary since entanglement cannot increase 
under SEP [8[. It has been shown that any deter- 
ministic controlled-unitary operator on two qubits 
implemented with bipartite LOCC using a resource 



of two entangled qubits necessarily requires a max- 
imally entangled resource Our paper takes a 
different approach to the problem, using SEP, and 
provides a proof applicable to general unitaries of 
arbitrary dimension. We show that if the resource 
has Schmidt rank equal to that of the unitary, the 
resource must be uniformly entangled in the sense 
that all its nonzero Schmidt coefficients are the same 
[Theorem 1(b)]. These same restrictions apply to 
LOCC, as it is a particular case of SEP. 

It is not hard to see that if the Schmidt rank of 
the resource is greater than the Schmidt rank of the 
unitary, then the resource need not be uniformly en- 
tangled (e.g. a larger rank resource that is majorized 
by a smaller rank maximally entangled state). We 
have found that it is in fact possible for such a larger 
rank resource to have less entanglement than would 
be required for a resource of Schmidt rank equal to 
that of the unitary. We have found examples of pro- 
tocols in both SEP and LOCC which deterministi- 
cally implement a controlled phase operation using 
less than one ebit of entanglement. In this case the 
unitary has Schmidt rank two and the resource has 
Schmidt rank three. Although the nonlocal unitary 
protocol given in [l^] can with certain probability 
consume less than one ebit of entanglement 0, we 
believe that ours is the first example of carrying out 
such a protocol deterministically using less than one 
ebit of entanglement. 

The remainder of this article is organized as fol- 
lows. Section [ll] sets up the problem of bipartite 
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^ Although the protocol given in lid is deterministic in the 
sense that it always succeeds in a finite number of steps, it is 
probabilistic in the amount of entanglement required. For 
any nontrivial unitary there is a chance that the protocol 
requires usage of the |i/'ci2) state, which has one ebit of 
entanglement. Thus, if the protocol only has access to a 
state with less than one ebit of entanglement there is a 
nonzero probability that the protocol cannot be carried out 
successfully. 
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deterministic implementations of unitary operators 
using SEP. Section Hill provides the requisite back- 
ground regarding map-state duality [Tl], [T^] and 
atemporal diagrams |13l4l5l |. Our main result is 
proved in Sec. IIVI using what we believe to be a 
new method based on the use of map-state dual- 
ity. In Sec. |V] we consider the case of a resource of 
larger Schmidt rank. There is a brief conclusion in 
Sec. I VII An appendix details the implementation of 
a controlled unitary using a qutrit resource state of 
less than one ebit of entanglement. 



II. NONLOCAL UNITARIES VIA 
SEPARABLE OPERATIONS 

We are interested in carrying out a bipartite uni- 
tary map U : "Ha ®'Hb-^ "Ha ® "^S' using as a 
resource an entangled state \ip) on two ancillary sys- 
tems Ha and Hb, by means of a separable operation 
{Ek Fk}, fc 1, 2, . . .. Here Ek:nA<»na^ Ha 
and Fk : Hb ®'Hb ^ V-b together form a product 
Kraus operator. For U to be unitary it is neces- 
sary that the dimensions of the Hilbert spaces satisfy 
d-AdB — dAd-B, but we do not require that (Ia = dA 
or dA = ds- The separable operation must satisfy 
the usual closure condition [l6| 

^{Ek®Fk)\Ek<»Fk)^lA®Ia®h®lB (1) 
k 

which is depicted in Fig.[TJa). 

In addition, for | $) any pure input state on Ha ® 
Hb, the outcome of the operation will be a pure 

state 

C/(|$) (<i>|)[/t = 

Y,{Ek ® Fk) (|$) ($1 ® W (M) {Ek ® Fk^ (2) 

k 

on Ha <8) V-B ■ Since the protocol is assumed to be 
deterministic, every term on the right side is pro- 
portional to the same pure state and it must be the 
case that 

{Ek ® Fk) = akU, (3) 

with Qffc some complex number. Note that both sides 
of ([3]) are operators acting on Ha (8) Hb', Fig. OJa) 
will help interpreting it correctly. 

The resource li/j) is assumed to have a Schmidt 
rank of D^, which means it can be written in the 
form 

W =Y.K\a,)®\h). (4) 

i=l 



for suitable orthonormal bases {|ai)} and {|5i)} of 
Ha and Hb, with Schmidt coefficients > for 
i < D^. 

Similarly, the bipartite operator U is assumed to 
have a Schmidt rank of Djj, meaning that it can be 
written in the form [TtI 

Du 

U = ^^JL^A^(S)B^, (5) 
i=l 

where {Ai\ and {Bi} are bases of the operator spaces 
C{Ha,Ha,) and C{Hb,Hb,), orthonormal under 
the Frobenius (Hilbert-Schmidt) inner product, and 
fii > for i < Djj. Equivalently, Djj is the mini- 
mum number of terms needed in order to write U in 
the form 'Y^Ci®'Di, without requiring Ci or T^i to be 
from an orthonormal basis. 



III. MAP-STATE DUALITY AND 
DIAGRAMS 

Map-state duality [TTI . [T2j plays a central role in 
the proof that will follow. This is a general con- 
cept that is sometimes referred to as reshaping or 
a partial transpose [llj and in a specific manifesta- 
tion is known as the Jamiolkowski or sometimes the 
Choi-Jamiolkowski isomorphism. States and maps 
are considered to both be tensors, and when a choice 
of orthonormal basis is fixed there is a natural linear 
relation between bras and kets (i.e. \i) O {i\ for all 
basis vectors 0. 

With this identification between bras and kets in 
place, the bipartite state jV') on the Hilbert space 
Ha ® Hb can be identified with the linear map V'' '■ 
Hb — > Ha obtained by turning kets into bras on the 
Hb space: 

IV') = Y. 1"') ® l^j) ^ ^' = II 

ij ij 

Similarly, the operators U, Ek, and F^, give rise to 
U' : Hb ®Hb ^ Ha ® Ha (by turning bras into 
kets on Ha and kets into bras on Hq), E^ : Ha 
Ha ® Ha (by turning bras into kets on Ha), and 
F^,^ -.Hb ® Hb — > Hb (by turning bras into kets on 



^ It is also possible to formulate map-state duality in a basis 
independent manner , however this is not necessary for 
the present work. 
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Hb and kets into bras on Hb), 

=Y. {Aj,B^\U\A,,B^)\A,,Aj){B,^,B^\, 

ij 77171 

E'k = ^{Aj\Ek\A,,a„,) \A,,Aj) (a„| , 

ijm 

Ff = {BjmB.,,b^) \b^} B,\ . (7) 



In the case of these three operators, map-map du- 
ality may be a more precise term, however we will 
use map-state duality to refer to any such partial 
transpose. The primed operator for Fk is denoted 
as in order to draw attention to the fact that 
its domain and range are swapped in comparison to 
El 

The equations introduced so far make use of six 
distinct Hilbert spaces and tensors of various rank. 
In such situations the underlying structure of equa- 
tions can be somewhat hidden when expressed us- 
ing Dirac notation. Abstract index notation is more 
transparent but can become unwieldy. For this rea- 
son we provide atemporal diagrams, similar to those 
found in , which should aid the reader in follow- 
ing the arguments in the text. 

Operators are designated by squares or rectangu- 
lar boxes. As a matter of style, the state ji/') and its 
corresponding operator ^' will be represented as a 
circle instead of a square. Lines between these boxes 
represent tensor contraction, and these lines are la- 
beled by the Hilbert spaces which they correspond 
to. Open lines on the left of a diagram represent the 
input to the total linear operator defined by the dia- 
gram, and open lines on the right represent outputs. 
Putting the inputs on the left means that operators 
are to be applied in a left-to-right manner, oppo- 
site to how algebraic equations are interpreted. As 
has been so far described, our diagrams are to be 
interpreted in exactly the same way as traditional 
quantum circuits as used for example in Nielsen and 
Chuang [l8|. The primary difference between our 
diagrams and traditional circuits is that in the lat- 
ter the horizontal direction is understood to repre- 
sent the passage of time whereas our diagrams make 
no reference to time. The presence of a summation 
symbol has the obvious meaning: the linear opera- 
tor depicted in the diagram denotes the terms of a 
series. The trace or partial trace operation is just 
a special case of tensor contraction and is denoted 
by joining the relevant spaces with a line. The iden- 
tity operator is represented by a line. With minor 
changes in style our diagrams are equivalent to the 
atemporal diagrams of [13| , and resemble other such 
schemes [l4l. Iisll. 
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FIG. 1. Atemporal diagrams, explained in Sec. IIIII (a) 
Closure condition, ([T}. (b) Apply and simplify using 
the adjoint of Fig. [^a) to get Q. (c) Multiply on the 
right by U to get (|1U|) . (d) Apply map-state duality to 
get (e) Restrict spaces to supports and ranges of 

operators to get ([12|. (f) Multiply by [/"^ (g) Trace 
over Hg to get (fT3|) . 



IV. ENTANGLEMENT REQUIREMENTS 

Our main result is the following: 

Theorem 1. Suppose that a unitary operator U is 
implemented deterministically by a separable opera- 
tion that makes use of the pure state entanglement 
resource {ip) [i.e. suppose that ^ and ^ hold]. 
Then 

(a) The Schmidt rank of {ip) is greater than 
or equal to the Schmidt rank Du of U . 

(h) If the Schmidt ranks are equal, Djj = D^, then 
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FIG. 2. (a) Deterministic unitary operation, ([3]). (b) 
Apply map-state duality to get ([8]). (c) Restrict spaces 
to supports and ranges of operators to get (|14|l . 



1-0) must be a uniformly (maximally) entangled state: 
all the nonzero Schmidt coefficients are the same. 

Proof of (a). Making use of map-state duality and 
the operators defined in ([6]) and ([7]), equation (jS]) 
[Fig. mja)] can be rewritten as [Fig. [H^b)] 



E'k^'Fl,^ 



akU'. 



(8) 



Since the rank of a product of linear operators is 
at most the smallest of the ranks of the individual 
operators, it follows that rank('(/;') > rank(C/'). The 
rank of an operator is equal to the number of its 
nonzero singular values. Since the Schmidt decom- 
positions (HI) and ([5]) are essentially singular value 
decompositions of if)' and [/', it is apparent that 
rank(-i/;') = and rank(C/') = Du and the inequal- 
ity becomes > Djj. Part (a) is proved. □ 

Proof of (h). Apply the closure condition ^ to (-01 
and use the adjoint of ([3]) to obtain 



Y,alU\Ek®Fk) = (01® /as. 



(9) 



as shown in Fig. [Ijb). Next, multiply both sides on 
the left by U to arrive at 



Y,o^l{Ek®Fk) = {i^\®U, 



(10) 



as shown in Fig. [Ijc). Making use of map-state du- 
ality gives [Fig. ^A)] 



Y,al{Ei®Fl^)^^'^®U'. 



(11) 



The map U' may in general have rank less than 
the dimension of <8) "H^ or 'Ha ® ^-a (which need 



not be equal to each other). In this case it will be 
useful to denote by 'H^ the subspace oi Hb T-ig 
which forms the support (or co-image or row space) 
of [/', the orthogonal complement of its kernel (null 
space) , and by T-Lj^ the subspace of T-La <E) Ha that 
forms the range (or image) of U'. Each of these 
subspaces has a dimension equal to Du, and U' is a 
nonsingular (invertible) linear map of onto , 
which we hereafter denote by U. In the same way 
one can introduce subspaces and 'Ha of ^6 and 
T-la which form the support and range of 0', and 
define to be the corresponding nonsingular map 
of rank from Hf^ to fia- Next, Ek is E'f, with 
its domain restricted to Ha, which can be strictly 
smaller than the support of i?^., and with its range 
restricted to T-l^, which could be smaller than the 
image of E'^. . Finally, F/^ is F^.^ regarded as a map 
from Hb ® Hb to Hb, but with domain and range 
restricted to Hg and respectively 0. 

The result of restricting ([TT|) to the subspaces just 
defined is 



(12) 



corresponding to Fig. [Ije). Multiplying on the left 
by and tracing over gives 



(13) 



see Fig. [TJf) and (g). Restricting (|8]) to subspaces 
results in 



Ek^pF^ 



akU, 



(14) 



see Fig.m^c). Here we have restricted the spaces over 
which matrix multiplications are being performed 
(Hi and Ha instead of Hb and Ha), however equal- 
ity is still maintained because the dimensions which 
have been eliminated correspond to the zero Schmidt 
coefficients of |0), which is to say the zero singular 
values of ip'. 

To complete the proof, make use of the assump- 
tion D^, ~ Djj. Then is also the rank of Ek and 
F^: all four operators in ((T4l) are full rank. Taking 
the inverse of both sides and inserting the result for 
tj^^ in leads to the result 

Y,Wk\H-^ ^ i^-^ ^ D^^\ (15) 

k 

where X]fc — 1 follows from (P), ^ and the 

normalization of 1-0) • With |0) in Schmidt form, -0 is 



It is significant that we define as Fj^ restricted to 
subspaces. In general it is not the case that Fj. is equal to 
restricted to subspaces. 
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diagonal, so ■0 = // \/D^. Therefore all the nonzero 
Schmidt coefficients of jV') are equal to 1 / yJD^ . 

□ 



V. LARGER RANK RESOURCE 

We have proved that a resource that is of the 
smallest viable Schmidt rank must be maximally en- 
tangled, but it is also possible to use a resource that 
is of higher Schmidt rank that is not maximally en- 
tangled. For one thing, if such a state meets an 
appropriate majorization criterion it can be deter- 
ministically transformed into a maximally entangled 
state . In this case the larger rank initial resource 
would have greater entanglement than would be re- 
quired if the smaller maximally entangled state had 
been used in the first place. There is however the 
possibility that some protocol could be devised to 
use a resource of larger Schmidt rank that has less 
entanglement than the maximally entangled state of 
smaller rank. 

In fact, we have numerically found examples 
of such constructions in both SEP and LOCC. 
One solution in SEP uses a resource state = 
VO^IOO) + V0.095(|ll) |22)) on two qutrits to 
implement the two qubit controlled unitary opera- 
tor U = diag{l, 1,1,6^"^} with = 2cos-i(35/36). 
We have verified this to be an exact solution using a 
computer algebra system. This resource constitutes 
less than one ebit of entanglement: the Von Neu- 
mann entropy is approximately 0.89 ebits. Since en- 
tropy cannot increase under SEP [H it is necessary 
for the resource that is consumed to have greater en- 
tanglement than the entangling capacity of the uni- 
tary being implemented. The entangling capacity of 
this unitary is shown in (20| to be approximately 0.23 
ebits. Since this is much less than the 0.89 ebits that 
we use, there remains the possibility that a different 
construction or an even larger rank resource could 
potentially lower the entanglement cost further. 

We also found an LOCC protocol which, though 
less efficient than the SEP construction just de- 
scribed, allows one to carry out a bipartite uni- 
tary deterministically using a resource with less than 
one ebit of entanglement. The resource in this 
case is = V0^|00) -I- %/(U(|ll) + |22)) and the 
unitary implemented is U — diag{l, 1, 1, e**} with 
— G.087r. The Von Neumann entropy of this re- 
source is approximately 0.92 ebits, and this is a four 
round protocol (Ahce, Bob, Alice, Bob). 

The constructions described above are instances of 
a more general continuous family of solutions that 
we have found, covering a range of controlled phase 
operations. As should be expected, a larger phase 
(j) requires a larger entanglement resource. In both 
the SEP and the LOCC case only certain classes of 



solutions were searched for, so it is possible that a 
more thorough search would provide more efficient 
protocols. The details of our SEP construction are 
presented in Appendix [X] Our LOCC construction 
consists of a long list of Kraus operators in numerical 
form, which is available upon request. 



VI. CONCLUSION 

We have shown that a unitary operator of Schmidt 
rank D implemented as a bipartite separable oper- 
ation requires an entanglement resource of Schmidt 
rank at least D. If the Schmidt rank of the resource 
is exactly equal to D, the resource must be uniformly 
(maximally) entangled with equal nonzero Schmidt 
coefficients. These restrictions apply also to LOCC, 
which is a subset of SEP. The proof uses map-state 
duality in a way which has not (so far as we know) 
been previously applied to problems of this type, so 
might have other interesting applications. 

Numerical results show that the amount of entan- 
glement required for the resource can be lowered by 
using a resource of Schmidt rank larger than D. A 
four round LOCC protocol has been found which 
uses a two-qutrit resource state with less than one 
ebit of entanglement to implement a bipartite con- 
trolled phase gate (albeit with a small phase). 

Although some large classes of unitaries are known 
to have implementations in LOCC using resources 
having the minimal Schmidt rank required by The- 
orem 1(a) [H-ll, Q, it is not known whether such 
minimal-rank implementations are possible for all 
unitaries. Given a unitary of Schmidt rank Djj it 
is always possible to find a collection of operators 
{Ek <^ Fk] such that © and ^ are satisfied with 
a resource of Schmidt rank = Djj. But it is 
not known if there is a separable operation satis- 
fying both (P) and ([3]). Consequently, it is possible 
that some unitaries may require a resource of greater 
rank than the lower bound given in Theorem 1(a). 
Even if such a minimal rank solution is always pos- 
sible in SEP, it still might not be possible in LOCC. 
This stands in contrast to the case of SLOCC where 
it is known that any unitary can be implemented 
using a state of Schmidt rank equal to that of the 
unitary Q. 
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Appendix A: Less than one ebit in SEP 

We performed a numerical search for solutions 
to ([Ij and with the resource and unitary tak- 
ing the forms 



IV-) = + v/(l - co)/2(|ll) + |22)), (Al) 

U^dmg{l,l,l,e'0}. (A2) 

In this case the unitary U is Schmidt rank 2 and the 
resource is Schmidt rank 3, so the spaces Ha and 
Hb are each 3 dimensional. In order to reduce the 
search space we looked for operators {Ek} and {Fk} 
of the form 



Ek = E^Sk and K = F^T^ 



(A3) 



where Sk ■ Ha Tic, Tk : Hb ^ He with He being 
a two dimensional space, and 



£/* — 



F, = 



1 




AA 





1 

1 

e'" 



AA 



BB 



(A4) 



(A5) 



It is possible to take advantage of the symmetry of 
the resource \tp) by searching for operator sets of the 
form 

{S-fci'M^A^"} and {TtL^M^'N'^} (A6) 
where l,m,n ^ {0, 1} and L, M, and N are defined 

by 



(A7) 



M = diag{l, 1,-1}, and N = diag{l, -1, 1}. There 
is no loss of generality in this assumption, since if 
{Tk}) gives a solution to ([T]) and ^ then so 
does ({-i|S'feL'M"W"},{TfeL'M"iV"}). This de- 
creases the number of independent operators (in- 
dexed by k) that need to be solved for, and in fact it 
turns out to be sufficient to consider only two values 
of fc. 




Initially we searched for solutions with 6 = tt/4 
and Co = 0.6 which, although representing more 
than one ebit of entanglement, is not majorized by 
a fully entangled resource of Schmidt rank 2. Once 
a solution was found, the parameters were variated 
until a value of cq was reached which represented a 
resource of less than one ebit of entanglement. Fur- 
ther constraints were added and variations made to 
simplify the solution and identify relations between 
the parameters. A family of solutions was found of 
the form (IA6I) with 



(A8) 
(A9) 
(AlO) 
(All) 



^1 = 

To- 
Ti = 

P = 



-1 1 -p 

—X — y * * 

—X -|- y * * 

(a; + ?;)e"'''/2 * * 



1 - s 



(A12) 



s = (1 - cos {9/2)) +y^{l + cos (9/2)) , 

(A13) 

where the parameters x, y, cq, and 9 must be solved 
for numerically. The asterisks in To and Ti repre- 
sent parameters that can be found using the relation 
Sk^'T^ = //4. 

A sequence of solutions for x, y, cq, and 9 were 
fed into an inverse symbolic calculator of our own 
design which uses a lookup table to convert float- 
ing point numbers into algebraic expressions. One 
of these solutions produced particularly simple alge- 
braic expressions: 



X = 9/5, 
Co = 0.81, 

9 = 2arccos(35/36). 



(A14) 
(A15) 
(A16) 
(A17) 



With this algebraic solution in hand, we used the 
computer algebra package Sage [2l| to verify that 
this indeed represented an exact (not just approx- 
imate to within floating point precision) solution 
to H]) and (121). 



[1] J. Eisert, K. Jacobs, P. Papadopoulos, and M. B. 

Plenio, Phys. Rev. A 62, 052317 (2000). 
[2] B. Reznik, Y. Aharonov, and B. Groisman, Phys. 

Rev. A 65, 032312 (2002). 
[3] L. Yu, R. B. Griffiths, and S. M. Cohen, Phys. Rev. 

A 81, 062315 (2010). 



[4] S. M. Cohen, Phys. Rev. A 81, 062316 (2010). 

[5] V. Gheorghiu, L. Yu, and S. M. Cohen, Phys. Rev. 

A 82, 022313 (2010). 
[6] N. B. Zhao and A. M. Wang, Phys. Rev. A 78, 

014305 (2008). 
[7] W. Dur, G. Vidal, and J. I. Cirac, Phys. Rev. Lett. 



7 



89, 057901 (2002). 
[8] V. Gheorghiu and R. B. Griffiths, Phys. Rev. A 78, 

020304 (2008). 
[9] A. Soeda, P. S. Turner, and M. Murao (2010), 
arXiv:1008.1128. 

[10] J. I. Cirac, W. Diir, B. Kraus, and M. Lewenstein, 
Phys. Rev. Lett. 86, 544 (2001). 

[11] K. Zyczkowski and I. Bengtsson, Open Syst. Inf. 
Dyn. 11, 3 (2004). 

[12] I. Bengtsson and K. Zyczkowski, Geometry of Quan- 
tum States (Cambridge University Press, Cam- 
bridge, 2006). 

[13] R. B. Griffiths, S. Wu, L. Yu, and S. M. Cohen, 
Phys. Rev. A 73, 052309 (2006). 

[14] G. E. Stedman, Diagram Techniques in Group The- 
ory (Cambridge University Press, 2009), 1st ed.. 



ISBN 9780521119702. 

[15] P. Cvitanovic, Group Theory: Birdtracks, Lie's, and 
Exceptional Groups (Princeton University Press, 
2008), ISBN 9780691118369. 

[16] R. Horodecki, P. Horodecki, M. Horodecki, and 
K. Horodecki, Rev. Mod. Phys. 81, 865 (2009). 

[17] J. E. Tyson, Journal of Physics A: Mathematical 
and General 36, 10101 (2003). 

[18] M. A. Nielsen and I. L. Chuang, Quantum Compu- 
tation and Quantum Information (Cambridge Uni- 
versity Press, 2000), 1st ed., ISBN 9780521635035. 

[19] M. A. Nielsen, Phys. Rev. Lett. 83, 436 (1999). 

[20] M. S. Leifer, L. Henderson, and N. Linden, Phys. 
Rev. A 67, 012306 (2003). 

[21] W. Stein et al.. Sage Mathematics Software (Ver- 
sion 4- 6-1), The Sage Development Team (2010), 
|http : / /www ■ sagemath . org 



